Nature of magnetic excitations in superconducting BaPe1.9Nio.1As2 
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Since the discovery of the metalhc antiferromagnetic (AF) ground state near superconductivity in 
iron-pnictide superconductors [1—3] , a central question has been whether magnetism in these materials 
arises from weakly correlated electrons [4, 5], as in the case of spin-density- wave in pure chromium [6], 
requires strong electron correlations [7], or can even be described in terms of localized electrons [8, 
9] such as the AF insulating state of copper oxides [10]. Here we use inelastic neutron scattering 
to determine the absolute intensity of the magnetic excitations throughout the Brillouin zone in 
electron-doped superconducting BaFe1.9Nio.1As2 (T^ = 20 K), which allows us to obtain the size of the 
fluctuating magnetic moment (m^), and its energy distribution [11, 12]. We find that superconducting 
BaFei 9Nio.iAs2 and AF BaFe2As2 [13] both have fluctuating magnetic moments (rn?) ss 3.2 /i^ per 
Fe(Ni), which are similar to those found in the AF insulating copper oxides [14, 15]. The common theme 
in both classes of high temperature superconductors is that magnetic excitations have partly localized 
character, thus showing the importance of strong correlations for high temperature superconductivity 
[16]. 

In the undoped state, iron pnictides such as BaFe2As2 form a metallic low-temperature orthorhombic phase with 

the AF structure as shown in Fig. la [17]. Inelastic neutron scattering measurements have mapped out spin waves 
throughout the BriouUion zone in the AF orthorhombic and paramagnetic tetragonal phases [13]. Upon Co- and 
Ni-doping to induce optimal superconductivity via electron doping, the orthorhombic structural distortion and static 
AF order in BaFc2As2 are suppressed and the system becomes tetragonal and paramagnetic at all temperatures [18]. 
In previous inelastic neutron scattering experiments on optimally electron-doped Ba(Fe,Co,Ni)2As2 superconductors 
[11, 12, 19-22], spin excitations up to ~120 meV were observed. However, the lack of spin excitations data at higher 
energies in absolute units precluded a comparison with spin waves in undoped BaFe2As2. Only the absolute intensity 
measurements in the entire Brillouin zone can reveal the effect of electron-doping on the overall spin excitations spectra 
and allow a direct comparison with the results in the AF insulating copper oxides [14, 15]. For the experiments, we 
chose to study well-characterized electron-doped BaFe1.gNio.1As2 [20, 22] because large single crystals were available 
[23] and their properties are similar to Co-doped BaFe2As2 [11, 12, 19, 21, 24]. 

By comparing spin excitations in BaFe1.9Nio.1As2 and BaFe2As2 throughout the Brillouin zone, we were able to 
probe how electron-doping and superconductivity affect the overall spin excitations spectra. We demonstrate that 
while the low-energy spin excitations are affected, the high-energy excitations show a very weak temperature and 
doping dependence. Comparison of our results with various theories suggests that neither a fully itinerant nor a 
localized picture explains the magnetic excitation spectrum. However, a combination of density functional theory 
(DFT) and dynamic mean field theory (DMFT) provides a natural way to improve on both these pictures. 

Figures Ic-e summarize our key findings for the electron-doped iron arsenide superconductor BaFe1.9Nio.1As2 and 
the comparison with the spin waves in BaFe2As2. The data points in Figs. Ic and Id show the dispersion of spin excita- 
tions for optimal doped BaFe1.9Nio.1As2 along [1, K] and [H, 0], and the solid lines show the fit of BaFe2As2 spin waves 
to an effective Heisenberg Jia — J\h — J2 model with Ji^ ^ Jif, [13] . Figure le shows the local dynamic susceptibility 
per formula unit (fu.), which contains two Fe(Ni) atoms, in absolute units, defined as = J x"(q, a;)(iq/ j dq 

[12], where x"{'i^^) = (l/3)i?'(Xa,3(l' '^))' different energies for BaFe2As2 and BaFe1.9Nio.1As2. It is clear that 
electron doping on BaFe2As2 only affects the low-energy spin excitations by broadening the spin waves below 80 meV, 
but has no impact on spin waves above 100 meV (see supplementary information). The quasiparticles that form within 
the spin-density- wave gap are sensitive to the Fermi surface change upon doping BaFe2As2, and hence the resulting 
low energy itinerant spin excitations substantially change, while the higher energy spin excitations are hardly affected. 

To substantiate the key conclusions from the data and calculations presented in Figure 1, we show in Figure 2 the 
two-dimensional constant-energy (E) images of spin excitations of BaFe1.9Nio.1As2 in the {H,K) scattering plane for 
several Brillouin zones at 5 K. In the undoped phase, spin waves in BaFe2As2 exhibits an anisotropy spin gap of 
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FIG. 1: Summary of neutron scattering and calculation results. Our experiments were carried out on the MERLIN time-of-flight 
chopper spectrometer at the Rutherford- Appleton Laboratory, UK [25]. We co-aligned 28 g of single crystals of BaFe1.gNio.1As2 
(with in-plane mosaic of 2.5° and out-of-plane mosaic of 4°). The incident beam energies were Ei = 20, 25, 30, 80, 250, 450, 600 
meV, and mostly with Ei parallel to the c-axis. To facilitate easy comparison with spin waves in BaFe2As2 [13], we defined 
the wave vector Q at (g^,, qy, q^) as {H,K,L) = {qxa/2Tr,qyb/2iT,qzc/2n) reciprocal lattice units (rlu) using the orthorhombic 
unit cell, where a — b = 5.564 A, and c — 12.77 A. The data are normalized to absolute units using a vanadium standard 
[13], which may have a systematic error up to 20% due to differences in neutron illumination of vanadium and sample, and 
time-of-flight instruments, (a) AF spin structure of BaFe2As2 with Fe spin ordering. The effective magnetic exchange couplings 
along different directions are depicted, (b) RPA and LDA+DMFT calculations of x"{'^) in absolute units for BaFe2As2 and 
BaFe1.9Nio.1As2. (c) The solid lines show spin wave dispersions of BaFe2As2 for Jia 7^ Jib along the [1, K] and [H, 0] directions 
obtained in Ref. [13]. The filled circles and upper triangles are spin excitation dispersions of BaFe1.9Nio.1As2 at 5 K and 
150 K, respectively, (d) The solid line shows low energy spin waves of BaFe2As2. The horizontal bars show the full-width- 
half-maximum of spin excitations in BaFe1.9Nio.1As2. (e) Energy dependence of x"{^) for BaFe2As2 (filled blue circles) and 
BaFe1.9Nio.1As2 below (filled red circles) and above (open red circles) Tc- The solid and dashed lines are guide to the line. 
The vertical error bars indicate the statistical errors of one standard deviation. The horizontal error bars in (e) indicate energy 
integration range. 



A = 9.8 meV [26]. On doping, the anisotropy spin gap disappears and spin excitations form transversely elongated 
ellipses that decrease in intensity with increasing energy [12, 21]. For energy transfers of £' = 10±3 (Fig. 2a), 33 ±3 
(Fig. 2b), 43 ± 3 (Fig. 2c), 60 ± 10 (Fig. 2d), and 81 ± 10 meV (Fig. 2e), spin excitations are peaked at the AF 
wave vector Q — (1,0) in the center of the Brillouin zone shown as dashed square boxes. As the energy increases to 
E = 113 ± 10 (Fig. 2f) and 135 ± 10 meV (Fig. 2g), spin excitations start to split along the K-direction and form a 
ring around the F point. Finally, spin excitations near the zone boundary at E ^ 157 ± 10 and 214 ± 10 meV form 
four blobs centered at Q — (1, 1). 

In order to determine the dispersion of spin excitations for BaFe1.9Nio.1As2, we cut through the two-dimensional 
images similar to Fig. 2 along the [1, K] and [H, 0] directions. Figures 3a-3f show constant-energy cuts along the 
[1, K] direction for = 25 ± 5, 55 ± 5, 95 ± 10, 125 ± 10, 150 ± 10, and 210 ± 10 meV. The scattering becomes 
dispersive for spin excitation energies above 95 meV. Figures 3g-3i show similar constant-energy cuts along the [H, 0] 
direction. The solid lines in the Figure show identical spin wave cuts for BaFe2As2 [13]. Since both measurements 
were taken in absolute units, we can compare the impact of electron-doping on the spin waves in BaFe2As2. At 
i5 = 25 ± 5 meV, spin excitations in superconducting BaFe1.9Nio.1As2 are considerably broader in momentum space 
and weaker in intensity than spin waves (Figs. 3a and 3g). Upon increasing the excitation energy to 55 ± 5 meV, the 



3 



3) 10±3meV b) 33±3meV C) 43 ± 3 meV 



2 

^ 1 
3 

— : 

V- 

-2 




-2-1 12-2-1 12-2-1 12 



^) 60±10meV ^) 81±10meV ^) 113±10meV 




-2-1 12-2-1 12-2-1 12 

H (r.l.u.) H (r.l.u.) H (r.l.u.) 



FIG. 2: Constant-energy slices through the magnetic excitations of BaFe1.9Nio.1As2 at different energies in several Brillouin 
zones. The images were obtained after subtracting the background integrated from 1.8 < H < 2.2 and —0.2 < K < 0.2. 
The color bars represent the vanadium normalized absolute spin excitation intensity in the units of mbarn/sr/meV/f.u and the 
dashed boxes indicate AF zone boundaries for a single FeAs layer. Two dimensional images of spin excitations at (a) E — 10±3, 
(b) 33 ± 3, (c) 43 ± 3, (d) 60 ± 10, (e) 81 ± 10, (f) 113 ± 10, (g) 135 ± 10, (h) 157 ± 10, and (i) 214 ± 10 meV. 



dispersive spin waves in BaFe2As2 become weaker and broader (Figs. 3b and 3h). For energies above 95 meV, spin 
excitations in BaFe1.9Nio.1As2 are almost indistinguishable from spin waves in BaFe2As2 in both the linewidths and 
intensity (Figs. 3c-3f, and 3i). Based on these constant-energy cuts, we show in Figs. Ic and Id the comparison of 
spin excitation dispersions of BaFe1.9Nio.1As2 (filled circles and horizontal bars) with those of spin waves in BaFe2As2 
(solid lines). Inspection of Figs. 1-3 reveals that electron-doping to BaFe2As2 only broadens and suppresses low 
energy spin excitations and has no influence for spin waves above 100 meV (see supplementary information). 

To reveal further the effect of electron doping on spin waves of BaFe2As2, we show in Figures 4a-4d constant-Q cuts 
at different wave vectors along the direction for spin excitations in BaFe1.9Nio.1As2. Near the Brillouin zone 

center at Q = (1, 0.05) and (1, 0.2), well-defined spin excitations are observed near E = 40, and 60 meV as shown in 
Figs. 4a and 4b, respectively. The intensity of the scattering from the spin excitations in BaFe1.9Nio.1As2 is, however, 
much lower than that of BaFe2As2 shown as solid lines in the Figures. On increasing the wave vectors to Q = (1, 0.35) 
and (1,0.5), the magnetic scattering peak near E = 100, and 120 meV, and are essentially indistinguishable from 
spin waves in BaFe2As2 as shown in Figs. 4c and 4d. Furthermore, spin excitations have virtually no temperature 
dependence between 5 K and 150 K (Fig. 4b). 

Finally, we show in Figs. 4e and 4f temperature dependence of spin excitations at energies near the neutron spin 
resonance E = 9 meV [20, 22] and at = 90±5 meV, respectively. While the intensity of the resonance at E = 9 meV 
increases dramatically below Tc, consistent with earlier work [20, 22], spin excitations at 90 ± 5 meV are identical on 
cooling from 150 K to 5 K. We note that high-energy spin waves in BaFe2As2 are also weakly temperature dependent 
[13]. Figure 4g shows the energy dependence of the dynamic spin-spin correlation lengths, which are about ^ « 14 
A and excitation energy independent. For comparison, the dynamic spin-spin correlation length (the solid line in 
Fig. 4g) in BaFe2As2 decreases with increasing energy and becomes similar to that of BaFe1.9Nio.1As2 for excitation 
energies above 100 meV. 

To check if spin excitations in the AF BaFe2As2 and superconducting BaFe1.9Nio.1As2 can be understood in an 
itinerant picture, we calculate the local susceptibility x"('^) using the random phase approximation (RPA) based 
on realistic Fermi surfaces and band structures [27]. Within RPA, the polarization bubble x° is computed from the 
density function theory (DFT) Kohn-Sham Green's functions while the irreducible vertex T"'' is approximated by the 
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FIG. 3: Constant-energy cuts of the spin excitation dispersion as a function of increasing energy along the [1, K] and [H, 0] 
directions for BaFe1.9Nio.1As2. The solid lines show identical cuts for spin waves of BaFe2As2 in absolute units, (a) Constant- 
energy cut along the [1, K] direction at £ = 25 ± 5, (b) 55 ± 5, (c) 95 ± 10, (d) 125 ± 10, (e) 150 ± 10, and (f) 210 ± 10 meV. 
(g) Constant-energy cut along the [H, 0] direction at -B = 25 ± 5, (h) 55 ± 5, and 95 ± 10 meV. The error bars indicate the 
statistical errors of one standard deviation. 



screened Coulomb parameters U and J. Using U — 1.3 eV and J = 0.4 eV and performing calculations above [27], 
we find that the RPA estimated x"{^) for BaFe2As2 and BaFe1.9Nio.1As2 (solid and dashed lines in Fig. lb) increases 
approximately linearly with energy and has absolute values about a factor of three smaller than the observation (Fig. 
le). Although the RPA calculation depends on Coulomb parameters used, we note that the 5-orbital Hubbard model 
calculation using U = 0.8 eV and J = 0.2 eV produces essentially similar local magnetic spectra [28]. Therefore, a pure 
RPA type itinerant model underestimates the absolute spectral weight of the magnetic excitations in iron pnictides. 

The solid blue and red lines in Fig. lb show the calculated local susceptibility using a combined DFT and DMFT 
in the paramagnetic state. Within DFT-fDMFT, x"(q, w) is computed by the Bethe-Salpeter equation using the 
polarization function x'^ and the two-particle local irreducible vertex function F*'"'' [27]. is computed from the 
interacting one-particle Green's function determined by the charge self-consistent full potential DFT-I-DMFT method 
and F""'' is extracted from the two particle vertex function of the auxiliary impurity problem. The latter is defined by 
the DMFT procedure using projection of all electronic states to the d character within the iron muffin-tin sphere. By 
comparing DFT-l-DMFT and RPA calculations in Fig. lb with data in Fig. le, we see that the former is much closer 
to the observation. Note that the calculation is done in the paramagnetic state, hence the low energy modifications 
of the spectra due to the long range AF order is not captured in this calculation. RPA can describe only the itinerant 
part of the electron spectra, while DFT-fDMFT captures the essential aspects of both the quasiparticles and the iron 
local moments formed by strong Hunds coupling (see supplementary information for more detailed discussion). The 
improved agreement of DFT-fDMFT thus suggest that both the quasiparticles and the local moment aspects of the 
iron electrons are needed to obtain the correct intensity and energy distribution of neutron scattering spectra [27]. 

One way to quantitatively compare spin excitations in iron pnictides with those in copper oxides is to estimate 
their total fluctuating moments, defined as (m^) = (Sh/n) J x" {^)du> / {1 — exp{—huj / kT)) [12]. Based on Fig. le, we 
find that (m^) = 3.17 ± 0.16 and 3.2 ± 0.16 fi% per Fe(Ni) for BaFe2As2 and BaFe1.9Nio.1As2, respectively. Using 
the formula for magnetic moment of a spin (m^) = {gfj,B)'^S{S -f 1) (where g = 2) [29], we find an effective iron spin 
S of about 1/2, similar to that of CaFe2As2 [30]. These results also show that superconductivity in electron-doped 
system hardly changes the total size of the fluctuating moment. In the fully localized (insulating) case, the formal 
Fe^"*" oxidation state in BaFe2As2 would give a 3dP electronic configuration. Hund's rules would yield S = 2 and 
(to^) — 24 /i^ per Fe. This is considerable more than the observed values suggesting that significant hybridization 
of Fe 3d with pnictide p orbitals and among themselves, which leads to a metallic state where the Hund's coupling 
is less important than in the atomic limit [31]. For comparison, we note that (to^) > 1.9 is^ per Cu for the AF 
insulating La2Cu04 measured over a similar energy range [14, 15]. From Fig. le, we see that the large fluctuating 
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FIG. 4: Constant- Energy/wave vector (Q) dependence of the spin excitations and dynamic spin-spin correlation lengths for 
BaFe1.9Nio.1As2 and BaFe2As2. (a)-(d) Constant-Q cuts at Q = (1,0.05), (1,0.2), (1,0.35), and (1,0.5), respectively, at T = 5 
(solid red circles) and 150 K (yellow filled circles) with background at Q = (2,0) subtracted. The negative scattering in the 
data are due to over subtraction of the phonon background. The solid lines are identical cuts from spin waves in BaFe2As2. For 
excitations below 100 meV, the intensity of the scattering of BaFe1.9Nio.1As2 is suppressed compared to that of BaFe2As2. For 
energies above 100 meV, the magnetic scattering is virtually identical between the parent and superconductor, (e) Constant- 
energy cuts at the neutron spin resonance energy of i? = 9 ± 1 meV [20] below and above T^. The solid lines are Gaussian 
fits on linear backgrounds, (f) Temperature dependence of spin excitations at £ = 90 ± 5 meV. (g) Energy dependence of 
the dynamic spin-spin correlation lengths (5) at 5 K obtained by Fourier transform of constant-energy cuts similar to those in 
Fig. 3a-f and Fig. 4e,f. For all excitation energies probed (10 < E < 200 meV), the dynamic spin-spin correlation lengths are 
independent of energy. The solid line shows energy dependence of ^ for BaFe2As2. The error bars indicate the statistical errors 
of one standard deviation. 



moment (m^) in iron pnictides arises mostly from high-energy spin excitations that is essentially temperature [13] and 
electron-doping independent within the errors of our measurements (Fig. 1). Since there are currently no high-energy 
spin excitation data in absolute units for optimally hole-doped Bao.67Ko.33Fe2As2 [32], it is unclear how hole-doping 
BaFe2As2 modifies the spin- wave spectra. 

The DFT-I-DMFT calculation suggests that both the band structure and the local moment aspects (e.g. Hunds 
coupling) of the iron electrons are needed to obtain a good description of the magnetic response in BaFe2As2 and 
BaFe1.9Nio.1As2. The weak electron-doping dependence of the fluctuating moment is consistent with the Hund's 
metal picture, where electron filling associated with the Fe 3dP electrons by Ni-doping is not expected to drastically 
affect the local moments. What is surprising is that the similarities in the local susceptibilities of the iron pnictides 
studied here and the parents of the cuprate superconductors. The large fluctuating moment, arising from Hund's 
rule coupling, and concentrated at higher energy in iron pnictides, nevertheless gives an imprint on the massive and 
anisotropic low-energy quasiparticles [33], which form Cooper pairs at low energy. This physics is different from the 
physics of doped charge transfer insulators appropriate for copper oxides [10], hence the electron correlations in iron 
pnictides and copper oxides have different microscopic origins, although they are important for understanding the 
magnetism and superconductivity for both materials. 
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Nature of magnetic excitations in superconducting BaFe1.9Nio.1As2 

Mengshu Liu, Leland W. Harriger, Huiqian Luo, Meng Wang, R. A. Ewings, T. Guidi, 
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Our BaFe2-xNia:As2 electron-doped samples were grown using self flux method as described before [1]. Since the 
electronic and superconducting properties of Co and Ni-doping of BaFe2As2 near optimal superconducting transition 
temperatures are almost identical [2], we chose to study spin excitations in optimally doped BaFe1.9Nio.1As2 with 
Tf, = 20 K. To further compare spin excitations in BaFe2As2 and BaFe1.9Nio.1As2, we show in Figure 1 spin excitations 
of these two materials at different energies. Consistent with data shown in Figs. 2-4 of the main text, we see that the 
effect of electron-doping is to mostly modify spin excitations below 80 meV. 

Our theoretical method for computing the magnetic excitation spectrum employs abinitio full potential 
DFT-I-DMFT method, as implemented in Ref. [3], which is based on the commercial DFT code of Wien2k [4]. 
The DMFT method requires solution of the generalized quantum impurity problem, which is here solved by the nu- 
merically exact continuous-time quantum Monte Carlo method (CTQMC) [5, 6]. The Coulomb interaction matrix for 
electrons on iron atom was determined by the self-consistent GW method in Ref. [7], giving {/ = 5eV and J = 0.8 eV 
for the local basis functions within the all electron approach employed in our DFT-I-DMFT method. 

The dynamical magnetic susceptibility x(q, w) is computed from the ah initio perspective by extracting the two- 
particle vertex functions of DFT-I-DMFT solution FJ^JJ. The polarization bubble x° is computed from the fully 
interacting one particle Greens function. The full susceptibility is computed from and the two-particle irreducible 
vertex function F^^JJ, which is assumed to be local in the same basis in which the DMFT self-energy is local, imple- 
mented here by the projector to the muffin-tin sphere [3]. In order to extract F^^^, we employ the Bethe-Salpeter 
equation (see Fig. 2) which relates the local two-particle Green's function {xioc), sampled by CTQMC, with both the 
local polarization function (x?oc) ^'^^ ^loc' 

Fj".' depends on three Matsubara frequencies (iv, ii/: iuj), and both the spin {(T\-i) and the orbital (ai_4) indices, 
which run over 3d states on the iron atom. T is the temperature. 

Once the irreducible vertex FJ^JJ is obtained, the momentum dependent two-particle Green's function is constructed 
again using the Bethe-Salpeter equation (Fig. 2) by replacing the local polarization function x?oc t>y the non-local one 

yO . : 

Finally, the dynamic magnetic susceptibility x(qj i^) is obtained by closing the two particle green's function with the 
magnetic moment /./, = fiB(L + 2S) vertex, and summing over all internal degrees of freedom, i.e., orbitals (ci!i_4), 
spins (cri_4) and frequencies {iv,iv'), on the four external legs 

X(q,iw)=Ty' V M°i-i/«°2-2 X"i<'i'''2-2 (iz/,w')q,iw (3) 

ivAv' "1"2 <'l<'2 
0'3^4 «^3«^4 

We note that the same abinitio methodology which is here used to compute the magnetic excitation spectra, was 
previously shown to describe the photoemission, the optical spectra and the magnetic moments of this material [8] in 
excellent agreement with experiment. 
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Figure 1: Constant-energy images of the spin excitations as a function of increasing energy for BaFei.gNio.i As2 and BaFe2As2 
in units of mbarns/sr/meV/f.u. (a) = 33 ± 3, (b) 43 ± 3, (c) 60 ± 10, (d) 81 ± 10, (e) 113 ± 10, (f) 135 ± 10, (g) 157 ± 10, 
and (h) 214 ± 10 meV. 
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Figure 2: The Feynman diagrams for the Bethe-Salpeter equation. It relates the two-particle Green's function (x) with the 
polarization (x") and the local irreducible vertex function (FJoc)- The non-local two-particle Green's function is obtained by 
replacing the local propagator by the non-local propagator. 
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